The Johnson-Segalman model with a diffusion term in cyhndrical Couette flow* 
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We study the Johnson-Segalman (JS) model as a paradigm for some complex fluids which are 
observed to phase separate, or "shear-band" in flow. We analyze the behavior of this model in 
cylindrical Couette flow and demonstrate the history dependence inherent in the local JS model. 
We add a simple gradient term to the stress dynamics and demonstrate how this term breaks the 
degeneracy of the local model and prescribes a much smaUer (discrete, rather than continuous) set 
of banded steady state solutions. We investigate some of the effects of the curvature of Couette 
flow on the observable steady state behavior and kinetics, and discuss some of the implications for 
metastabihty. 



1. INTRODUCTION 




ity, with cither a discontinuity or a discontinuous first 
derivative in the apparent constitutive curve. In some 
worm-like micellar solutions this flow behavior coincides 
with a range of strain r ates for which macroscopic phase 
separation is observed ( Berret et al. , 1994; pallaghan et 



1996 ; Grand et al. , 1997 ) under controlled strain rate 



applied strain rate y 

FI G. 1 . Local constitutive relation for the JS model 
(Eg. 3.1 ). The horizontal line shows a possible stress for 



conditions, with stresses on a well defined, reproducible, 
and history-independent plateau. 

Although wall slip or internal fracture, rather than 
bulk instability, has been suggested to explain the spurt 
effect (Dcnn, 199C), bulk phase separation is undisputed 
in the wormlike micelle system. This phenomenon has 
often been interpreted in terms of a non-monotonic con- 
stitutive curve analogous to a van der Waals loop, as 
in Fig. |l|. Here axy is the shear component of the 
fluid stress tensor, x is the coordinate in the flow di- 
rection, and y is in the gradient direction. There are 



many kinds of non-monotonic behavior (Porte, Berret 



banded flow, with material on either the low or high strain and Harden, 1997; Olmsted and Lu, 1997), and we focus 



rate branches. The negative slope portion of the curve repre- 
sents mechanically unstable homogeneous flows. 

Non-analytic flow behavior is seen in a wide range of 
complex fluids, i n cluding polymer m elts (Bagley, Cabot 
and West , 195S ; Vinogradovl, 1972), wormlike micelles 



here on systems with multiple strain rates for a given 
stress. Constitutive curves of this character are found 



( Johnson and Segalman 




1977 




Malkus, Nohel and Plohi , 


199C; 


Malkus, Nohel and Plohr, 1991), the Doi-Edwards 



and Smoulders , 1996; prand, Arrault and Cates , 1997 



Boltenhagen, Hu, Matthys and Pine , 1997 b and lamellar 



crrct, Roux and Port4 1994; pallaghan, Cates, Rofe wards , 1989| ) and the similar model of Cates ( 1990D for 



surfactant phases ( Roux, Nallet and Diat 
Roux and Nalletl |1993|; piat, Roux and 




Bierro and Roux "|1 997). In the "spurt e ffect" in pipe flow 



( Bagley et al. , 195S; Vinogradov , 1972) the flow rate in 



creases discontinuously above a given applied pressure 
gradient. In Couette flow this behavior manifests it- 
self in non-analyticities in the steady state flow curves 
of wall stress as a function of relative cylinder veloc- 



wormlike micelles (both augmented with a stabilizing 
high shear-rate bra nch), and other recent mo dels for 
branched polymers (McLeish and Larson, 199S). These 
flow curves have regions of negative slope, da^y/d'y < 0, 
which are hydrodynamically unstable. For applied strain 
rates in these unstable regions, a reasonable scenario con- 
sistent with experiments is for the system to break up into 
bands of different flow rates, each band lying on a stable 
part {daxy/dj > 0) of the curve. This possibility has 
been studied for various models of non- Newtonian fluids 
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( iMcLeish and Ball|, |1986|; iMcLeishj |1987|; l Olmsted and 



Gold lbarll, |1992|; [Spenley, Yuan and Gates] , |1996| ; plm 



sted l and Lu j_ 1997 ). However, the stress can apparently 
lie a nywhere within a range of stresses. This degeneracy 
of selected stresses has engendered much discussion, and 
is at odds with the wormlike micelles, which shear-band 
at a well-defined s t ress that is independent of flow his- 
tory ( iBerret et al.|, |1994 [Callaghan et al.| , [l996| ; [Grand 
et af], |1997| ; [Boltenhageii et al.| , |1997| )"! 

Most experiments are performed under conditions 
which impose inhomogeneous stress. For example, in 
pressure-driven (Poiseuille) pipe flow the steady state 



tional term to the JS model, which could be thought 
of as stress relaxation by diffusion of differently-strained 
polymer str a nds ( El-Kareh and Leal, 1989; Goveas and 



shear stress varies linearly with distance from the center 



of the pipe, while in cylindrical Gouette flow the stress 
vari es as the inverse of the square of the radius. There is a 



Fredrickson , 1998). Such "gradient terms" are not usu 
ally included in constitutive equations, but arise natu- 
rally in treatments of hydrodynamics of liquid crystals, 
binary fluids, and polymer blends, and are expected to be 
present (albeit small) even in one component melts. Such 
terms will be most important exactly in the highly inho- 
mogeneous interface between shear bands, and consti- 
tute a singular perturbation on the equatio ns of motion, 
with a profound effect on stress selectio n (p!^u, Olmsted 
and Ball, 199£; Radulescu and Olmsted, 1999[ ). Previous 
workers have used gradient terms to study stress selection 



gro'vdng body of theoretical work addressing shear band- and Goldbart 



in models for liquid c r ystals in planar sh e ar flo w ( Olmsted 



1992 ; Olmsted and Lu, 1997), and in a 



ing in inhomogeneous flow conditions; in particular, we 
are interested in the relationship between the degeneracy 
of models such as the local JS model, and the possible 
history-dependence in different flow geometries, and how 
to resolve this degeneracy. 

iVlcLeisn (1987) studied the stability and hysteresis of 



shear banding in a modifled local Doi-Edwards fluid in 
pipe flow, and found that normal st ress effects should af - 
fect the stability of bande d flows. Malkus et al, (199C) 



and 



Malkus et al 



(1991) studied the local JS model 
in planar Poiseuille flow and presented a comprehen- 
sive phase-plane analysis of the steady state behavior, 
showing how a hysteresis loo p resul ted from a particular 
flow history. Greco and Ball (1997) studied shear band- 
ing in start-up flows of the JS model in Gouette flow, 
and demonstrated that the curved geometry of Gouette 
flow produces a simple well-controlled banded flow, in 
contrast to planar flow, in which start-up calculations 



of s: milar local models give an uncontrolled number of 



bands. More recently the degener acy of the JS mode 



was stud i ed in cylindrical Gouet te ( Gcorgiou and Vlas- 



sopoulos , 1998 ) and Poiseuille (Fyrillas, Gcorgiou and 



VlassopouloE , 1999) geometries, with the conclusion in 
each case that selection of a given banded steady state 
depends on the initial perturbation. 

In this paper we first study the consequences of degen- 
eracy for the history dependence of banded solutions to 
the local JS model in cylindrical Gouette flow by impos- 
ing several flow histories. We flnd, similar to Malkus et 



al. (1991) and Fyrillas et al. (1999) with the JS model 



in pipe flow, that the apparent flow curves (in this case, 
wall stress as a function of relative cylinder, or "gap", 
velocity) depend on flow history. We restrict ourselves to 
four basic flow histories and demonstrate how to produce 
a continuum of possible solutions. Even in the curved ge- 
ometry, the solution is unstable with respect to noise in 
the initial conditions, a s in the planar case (Gcorgiou 
and Vlassopoulos , 1998 ). This is in distinct contrast to 
the many experiments which have found a reproducible 
and unique selected shear stress in systems as diverse as 
worm- like micelles and surfactant lamellar phases. 

To explore unique stress selection we add an addi- 



one-dimensional scalar model for banded flow in micelle 



solutions by Spenley et al. (1996). Although we make no 
claims that the JS model accurately approximates these 
systems, it is a simple model which possesses the requisite 
non-monotonic constitutive flow curve and whose dynam- 
ics, as suggested by Berret, is compatible to some aspects 
of kinetics experiments of shear banding in solutions of 
wormlike micelles ( p3erret et al. , 1994| ; Berret, 1997). 

The outline of this work is as follows. In Section ^ 
we present the equations for the JS model in Gouette 
flow, including the new diffusion term. In Section ^ we 
present the results from numerical experiments for pre- 
scribed strain rates (gap velocity) and different flow his- 
tories in the absence of diffusion; and in Section ^ we 
examine the effect of diffusion, which renders the model 
more realistic with respect to the body of work on com- 
plex fluids which display unique stress selection. We fin- 



ish with a discussion. In a companion paper (Radulescu 



and Olmsted, 1999) we analyze the diffusive JS model 
in Poiseuille and Gouette geometries, using asymptotic 
matching techniques, to analytically support many of the 
findings in this work. 



2. MODEL 



2.1. Dynamical Equations 



The Johnson-Segalman model has been discussed by 



many authors (lohnson and Segalman, 1977 


Malkus et 


al. 


, 1990; 


Malkus et al., 1991; Greco and Ball, 


1997). The 



momentum balance condition for an incompressible fluid 
is 



Pidt 



v-V) V = V-T. 



(2.1) 



where p is the fluid density, v the velocity flcld, and T 
the stress tensor. The stress tensor is given by 



T = -pI + 2r;D + S, 



(2.2) 



where the pressure p is determined by incompressibility 
(V-v = 0), rj is the "solvent" viscosity, and S is the 
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"polymer" stress. We separate the velocity gradient ten- 
sor (Vv)q^ = daVfs into symmetric and anti-symmetric 
parts D and ft, respectively: 



2D = Vv + ( Vv)^ 
2r2 = Vv - ( Vv)^ 



(2.3a) 
(2.3b) 



The non-Newtoni an "polymer" stress is ta k en to have the 
form proposed by Johnson and Scgalman (1977): 



♦ 



2^D-is, 

T r 



(2.4) 



where fi is the "polymer" viscosity, r is a relaxation time, 
and V is the diffusion coefficient. The difference from the 
usual JS model is the added diffusion term. This has been 



derive d due to dumbbell diffusion by El-Karch and Leal 
(1989), and is also expected to have non-local contribu- 
tions due to the spatial extent of macromolecules (Lu ct 
al. , |1999| ). Gradient contributions to the polymer stress 
of blends due to the latter effect have been derived by 
poveas and Fredricksoii| ( 1998| ). The time evolution of 
S is governed by the Gordon-Schowalter time derivative 
( pordon and Schowalter] , |1972| ), 

S = (Sf + V- V) S - (JIS - S J7) - a (DS + ED) . (2.5) 

The "slip parameter" a was interpreted by Johnson and 
Segalman as a measure of the non-affinity of the polymer 
deformation; i. e. the fractional stretch of the polymeric 
material with respect to the stretch of the flow field. For 
\a\ < 1 the polymer "slips" and the steady-state flow 
curve in planar shear is capable of the non-monotonic 
behavior in Fig. |l]. 

In this work we consider non-inertial flows (zero 
Reynolds number limit), for which the momentum bal- 
ance Eq. (2.1) becomes V-T = 0. Although inertial effects 
should not affect stationary solutions, they may influ- 
ence transients at the very early stages of start-u p exper- 
iments, for which inertia cannot be neglected. Spenley 
et a[] ( |1996| ) showed that inertial terms have practically 
no influence on the later stages of the dynamics, so at 
timescales of order the characteristic relaxation time r 
the transients are co rrectly described by th e non-inertial 
dynamics. Following Greco and Ball ( 1997|) , we assume a 
flow field with azimuthal symmetry, v = v{r,t) between 
concentric cylinders of radii Ri < i?2, in cylindrical coor- 
dinates {r, 9, z}. Integrating the balance condition yields 



r]j{r,t) + T.re{r,t) 



where we define the local strain rate by 



d_ 

dr 



(2.6) 
(2.7) 

(2.8) 



the integration constant F is the torque per cylinder 
length, and a{r,t) is the inhomogeneous shear stress in- 
side the Couette cell gap. We consider no-slip boundary 
conditions on the fluid velocity v, a fixed outer cylinder, 
and an inner cylinder moving at a prescribed gap velocity 
V. This implies a global constraint 



that may be used to find the torque. 



2 o2 



Ri 



dr 



(2.9) 



(2.10) 



-Ri 



2.2. Rescaled equations 



In cylindrical coordinates the dynamical equations for 
S, Eq. (H^J), become 



22? 

C T^rr = -(1 - a)'jY.re + (See - '^rr) 

2V 



Cllee = (l + a)7Ere 5- (Dee - Srr) 



/•V 41? 

LY^re = -7 T^re - 7 



1 -a. 



1 a. 



where 



(2.11a) 
(2.11b) 
(2.11c) 



(2.12) 



Next we change variab les. We introdu c e the space vari- 
able X G (0, 1) used by preco and Bali| ( |1997| ), 



1 ^2 
q = In — , 

^ Ri 



and deflne 



Z = 

W 
S 



l_a, 
2 

1 - a 



1 + a. 



2 



y ^ + V 



(2.13a) 
(2.13b) 

(2.14a) 

(2.14b) 
(2.14c) 



We now introduce the following dimensionless variables, 
where X is the dimensionless version of X: 



V = V 

z 



qRi 

T 



w = w- 
v = v 



q^Rl 



^iRl 



(2.15a) 
(2.15b) 
(2.15c) 
(2.15d) 
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After rescaling, Eqs. ( ^.11 ) become 



{W I aZ)a 



with 



C,W = ^ —{W + aZ) 

1 — 



(2.16a) 
(2.16b) 
(2.16c) 

(2.17) 



The local strain rate and torque (Eqs. 2^ 2.1C| ) may 
be written as 



= f - S 
2(7 



1 



-2q 



((5)-.V-), 



where 



(5) = dxS 
Jo 



1 



dr 



7- 



and 



(2.18) 
(2.19) 

(2.20) 
(2.21) 

(2.22) 



is the viscosity ratio. 

We consider the dynamics of flows for prescribed 
histories of inner cyhnder velocities V, by solving 
Eqs. ( |2.16y2.18| , |2l9| ). The only parameters are V, which 
sets the length scale of any interfaces; the slip parameter 
a; the viscosity ratio e; and the curvature of the Couette 
cylinders, determined by q. We will often use the relative 
gap size, 



p = e"? - 1 

i?2 — ^1 

" Ri 



(2.23) 
(2.24) 



instead of q. 



For the interpretation of the stress diffusion ter m as aris- 
ing from the di ffusion of polymeric dumbbells ( El-Kareh 



and Lea], 1989), this corresponds to a zero flux boundary 
condition. In a more realistic model gradient terms are 
expected to arise from inter- and intra-molecular interac- 
tions as well as center-of-mass diffusion, and the bound- 
ary condition may be a fixed non-zero gradient, or per- 
haps a fixed value for the stress. Clearly the boundary 
conditions depend on the detailed physics. In this work 
the nature of the boundary condition is expected to play 
a role in the limit when the interface touches the wall, so 
predictions in this limit may be non-universal. 



COUETTE FLOW WITH NO DIFFUSION: 

v = o 



3.1. Steady State Solutions 



The steady states and stability of the local JS model 
{V = 0) in planar shear flow have been studied by 
many authors. Although some numeri cal calculations 
have shown well-defined shear banding (Espahol, Yuan 
Bal] 



and 



1996; Yuan, 1999), it is generally believed 



that the local model supports s hear ban d s at any stress 
in the non-monotonic region (Renardy , 1995 ; Spenley 



al 



ct al,^1996; Georgiou and Vlassopoulo^ , |1998| ; |Lu et 
1999D . he presence of the diffusion term stress 



is selected at a well defined value, independent of his- 
tory. We will not discuss the case of 2 ? 7^ in p l anar fl ow 

(|l99|). 



in detail here, but refer the reader to Lu et al 



Greco and Ball ( |1997| ) performed start-up calcula- 
tions for I? = in a cylindrical Couette ge ometry and 
found, unlike the 1+lD planar simulations ( ^penley et 



al., ^996; Espafiol et al., 1996), an apparently repro- 



ducible steady state with two shear bands and a single 
zero width interface between bands. Before studying this 
model in detail, we present a graphical construction of 
the steady state banded solutions. 

The steady state solutions to Eqs. ( ^.16| ) for P = 
have the lo cal form as the flat case, given by (Greco and 



BaU 



19971) 



2.3. Boundary Conditions 

While the fluid velocity it taken to obey no-slip bound- 
ary conditions, for S we choose the following boundary 
condition: 



= 0, 



(2.25) 



Z = Q 
w ==^S 

S ■ 



7 



1 



■7 



(3.1a) 
(3.1b) 

(3.1c) 



Using the local force balance, Eq. (p.6D, and Eq. |3T| c) 
we obtain: 



d^w = d^z = d:,s = 0. 



(2.26) 



f^7 



1 + 7 



(3.2) 
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where three-dimensional space spanned by {F, 7, r} this relation 

- _ Ci 'x\ defines a surface, shown in Fig. [^a. Imposing the uniform 

Ri \ ■ I torque condition (the intersection of the surface r(r, 7) 

Eq. (IJ) defines the locus of steady states parametrized ^i*^!^ ^ plane) yields a local relation between strain-rate 

by torque, radial position, and local strain rate. In the and radius 7(r,t), shown in Fig. |b,c. 



(a) 




« 

c 

'cs 




(C) 



InR, 



In r„ In r. 



In R, In R, In 



Inr^ 



In R, 



FIG. 2. (a) The curved surface sho ws t he torque T{r, 7) as a function of strain-rate 7 and radius r for the Johnson-Segalman 
model in Couette flow, given by Eq. (3.2). The plane is at a specified torque. (b,c) Banding profiles at different torques for the 
JS model in Couette flow, with no diffusion term [V = Q), obtained from the intersection of the plane and curved surfaces in 
(a) according to Eq. ([3.2|). Thin lines are strain-rate 7 as a function of radius r; thick lines are banding solutions, (b) Possible 
profiles for torques Fa and Tb > ^ a at the same gap velocities V , given by the area of the shaded region, (c) Two possible 
bands for a given torque, again with the same gap velocity. 



A steady state banded flow profile traces a trajectory 
7(r) in Fig. with a vertical jump. In the two solu- 
tions shown in Fig. |^b, the high strain rate phase lies 
nearest the inner cylinder. Because of the inhomogeneity 
of Couette flow the high and low strain rate phases are 
non-uniform. 

Given a solution, e.g. the solid line in Fig. ||b for 
torque F^, t he shaded area corresponds to V , according 
to Eq. ( ^.2l| ). By placing the interface at different posi- 
tions, we can generate a continuum of solutions with dif- 
ferent gap velocities for a given torque. Conversely, it is 
easy to find a continuum of torques for a given gap veloc- 
ity, with each torque corresponding to a unique interface 
position. Upon increasing the torque to F^ > Fa while 
keeping the interface fixed (specified by the vertical solid 
line), the area under the new 7(r) curve, and hence the 
gap velocity V ^ increases. To recover the V we must move 
the interface closer to the inner cylinder (i?i), obtaining 
the solution shown in the thick dashed line. Also, mul- 
tiple interface solutions are possible: flow proflles which 
traverse between high and low strain rate branches an 
odd number of times are possible steady states. Hence, 
for a band sequence the steady state is completely defined 
by any two of the three quantities F, V/Ri^ and inter- 
face position r» . If both the low and high strain branches 
span the gap (Fig. ^), which happens for a small cur- 
vature {R1/R2 close to 1), both the conventional band 
sequence, and an inverted band sequence, with the high 
strain rate material at the outer cylinder, are possible. 
This possibility is stabilized for flatter cylinders, and in 
the planar limit the conventional and inverted bands are 



symmetry-related partners. 

Hence, a prescribed gap velocity does not select a 
unique stress (torque) in cylindrical Couette geometry. 
At flxed V, the family of single interface steady flow so- 
lutions is parametrized by F or by G [i?i,i?2]. In 
principle, a universal banded steady flow might still ex- 
ist if all the initial conditions or flow histories evolved to 
the same steady state configuration (selected values of r» 
and F). To explore this we next simulate flows prepared 
by different flow histories. 



3.2. Dynamics 



1. History dependence 



We have evolved Eqs. ( ^.16 ) through various sequences 
of gap velocities {V{ti)}, where ti is the time when the 
gap velocity is changed, and considered four types of flow 
histories: 

1. Jump-up: jump each time from a gap ve- 
locity Vdown, small and on the low strain 
rate branch, to a gap velocity V^; {V"(0) = 

Vdown, Vih) = Vi,V{t2) - Vdown,- ■■ ,V{t2n) = 
Vdown, V{t2n+l) = Vi}. In particular Vdown = 

means starting from rest. 

2. Jump-down: jump each time from a velocity Vup, 
large and always on the high strain rate branch, 
to a velocity V; {V{0) = Vup, V{ti) = Vi, T>(t2) = 



V 



up, ■ 



,Vit2n)^V^p,Vit2n+l)=Vn}. 
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3. Ramp-up: increasing sequence of velocity values; 
{l7(0) = < Vih) = < . . . < V{t^) = 14}. 

4. Ramp-down: decreasing sequence of velocity val- 
ues; {1/(0) = K > V{h) = Vn-l > ... > 

The lengths U — U-i of the relaxation intervals were 
chosen to be much larger than the characteristic re- 



laxation time for each velocity change, in order to ob- 
tain the steady state at the end of each interval. We 
use an implicit Crank-Nicholson algorithm to perform 
the dynamics. At each time step we must determine 
the torqu e T through the non-local int egral condition, 
Eq. (^.19 ). For comparison, we note that Greco and Ball 
(1997) performed a similar calculation using e — 0.05 and 
p = 0.1,0.02,0.01, and starting from rest. 
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FIG. 3. Torque vs. gap velocity for the four scenarios described in the text, for "D = and e = 0.05. a.) p — 0.1; b) 
p — 0.01. Vdown = 0.0, Vup = 12.0, the spatial mesh is 200, the time step is 0.005, and 10000 time steps assured steady state at 
each value for V . Ramps are shown for both large [dV > 0.2) and quasi-continuous {dV ~ 0.01, denoted qc) velocity changes. 



The steady state values of f at a given gap velocity 
V are history dependent, as shown in Fig. ^ For the 
same gap velocity, steady flows obtained via different 
scenarios differ by the number of bands and the pro- 
portion of high shear rate band. For high curvatures 
steady banded flow has two bands for all scenarios, with 
the high shear rate band at the inner cylinder, but the 
amount of high shear rate band for the same gap veloc- 
ity is maximum (the torque is minimum) on ramp-down, 
minimum on ramp-up, and takes different intermediate 
values for jump-up and jump-down scenarios. The case 
of low curvature p — 0.01 provides a surprise in jump- 
up from rest, also noted by Greco and Ball. Within the 
interval 1.96 < y < 5.1 (between points Ml and M2 in 
Fig. ^d) the steady state contains two interfaces, sepa- 
rating an interior high strain rate band from two outer 
low strain rate bands against the cylinder walls. For 
0.95 <V < 1.96 (between points B and Ml in Fig.|b) 
we obtain the usual high-low two-band sequence, while 
for 5.1 < f < 6.48 (between points M2 and O in Fig.|b) 
the two-band sequence is inverted (low-high) with the 
low strain rate band of material at the inner cylinder. 
The three-band region is separated from the two-band 
region by discontinuities of the slope, dT /dV (points Ml, 



M2). The values of the gap velocity separating these dif- 
ferent types of steady flow correspond to local minima 
of the torque T{V) (Fig |b). The time development of 
two-band and three-band profiles can be seen in Fig. IJ. 
The other three scenarios produce the normal two band 
sequence high-low. 

Finally, another history-dependent feature occurs on 
ramp-up and ramp-down scenarios. The size of the high 
shear rate band in steady banded flow and therefore the 
torque value depends on the magnitude of the differences 
y[ii+i) ~ between succesive velocities. If the ve- 

locity is changed quasi-statically (infinitesimal changes 
of V), flow curves are smooth and correspond to top and 
bottom jump (the total shear stress at the position of 
the interface corresponds to the maximum and minimum 
values of the local constitutive relation) on ramp-up and 
ramp-down respectively. Large velocity changes produce 
spikes in the flow curve, both on ramp-up and ramp- 
down. The reason for this is the pinning of the interface 
at a stress value which is lower than the maximum al- 
lowed value or higher than t he minimum al l owed value, 
and has been also noticed by Spenley et al. ( 1996 ) for a 
different diffusionless model. 
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2. Noise effects 

The dynamical selection of steady states in the four 
scenarios is unstable with respect to noise in the initial 
conditions. To show this we used perturbed initial con- 
ditions at each new velocity value, in the form 



Ep'(r) = s(i\(r;^)+ry(^)(r), 



(3.4) 



and breaks up into many bands whose widths are lim- 
ited only by the mesh size. The size of the unstable re- 
gion increases with the amplitude "qmax of the noise and 
it is smaller for higher curvature {p = 0.1). Diffusion- 
less dynamics of the banded flow is mor e stable in high 
curvat ure Couette cell, as suggested by [Greco and Ball 
(|T997|). 



where (r; oo), j = 1,3 are the three steady- 

state components of the material stress obtained at 
step i-1, and rj^^\r) are random, uniformly dis- 
tributed in the interval {—rjmaxiVmax) and uncorrelated, 
lri^^)(r)ri^^'){r')) = for r 7^ r', j 7^ /. Fig. || shows the 
result of two-state experiments (z = 1,2), starting from 
rest (state i = 1). A region in the high strain rate band 
close to the position of the interface (that would be se- 
lected without noise) is unstable with respect to noise. 



3. Basins of attraction and continuous degeneracy of steady 
handed flows 



The dynamical equations ( ^.16 ) with V = correspond 
to an infinite dimensional nonlinear dynamical system. 
The dynamics in the functional space of banded flows has 
a complex structure of attractors (steady banded flows). 
Without noise, certain classes of trajectories, correspond- 
ing to initial conditions inside basins of attraction, reach 
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the same steady state. The jump-up branch BG (Fig. ||) 
is a set of attractors, collecting trajectories starting on 
the low shear rate branch. For an imposed gap velocity 
Vend in the banded regime, jumping up from an initial 
homogeneous steady state of gap velocity Vdoum < Vp 
(initial state below the point F on the flow curve, Fig. ^a) 
leads to an identical final banded steady state (point Ei ) . 
For higher initial gap velocities Vp < Vdown < Vb (initial 
states between points F and B on the flow curve. Fig. ^a) 



the final torque spans the interval Ei — E2 between the 
jump-up and jump-down branches (the discreteness in 
Fig. H is a finite mesh effect). The jump-up from rest 
branch is thus the lower limit of a continuous family of 
jump-up attractors. Similarly, jumping-down from a ve- 
locity Vup > Vg ends on the attractor branch BC which 
represents the upper limit of jump-down attractors which 
collect trajectories with Vup < Vq (Fig. ||). 



1.2 




Time 



1.0 




Despite the interesting structure of these dynamics, ex- 
periments on wormlike micelles invariably display unique 
ban ding flows, independent of history and initial condi- 



tion 



This is incompatible with the numerical results 



obtained for I? = 0; as we shall see in the following, 



the imposition of non-vanishing stress diffusion I? ^ 
provides a route for constructing unique banding flows 
(lEl-Kareh and Lca| |198g| ; [Olmsted and Lu| , |1997] ; |Lu et 



al. 
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FIG. 6. Structure of basins of attraction. The arrows lines connect start and end states (attractors) upon (a) Jumping up 
from the low shear rate band (Vend = 1.75) (b) Jumping down from the high shear rate band (\4nd = 5). 



4. COUETTE FLOW WITH DIFFUSION: O / 



4.1. Flow curves and History independence 



V=3.0, p=0.1,E=0.05 
a=0.5, D=0.001 




the interface, with a selected stress Gsei independent of 
imposed gap velocity, interface position, curvature, or T). 

In contrast to planar flow, for which the flow curve 
has a flat stress plateau, in cylindrical Couette flow the 
torque in the banding part of the flow curve has a slight 
dependence on the imposed gap velocity (Fig. H). This 

is a consequence of the relation F — Ogei ft , with asei 
the selected stress at the interface position f, , the lat- 
ter which increases from 1 to 1 + p as the fast-flowing 
band gradually fills the gap with increasing gap velocity 
V. The slope of the plateau, dT/dV ~ p, vanishes in 
the planar limit of a thin gap. This provides a second 
mechanism for a sloped plateau, in addition to the ef- 
fects of a co ncentration difference between t he banding 



y constitutive curve 



sted and Lu 



states T{V) ( jchmitt, Marques and Lequeux , 1995 ; 01m- 



1997) 




FIG. 7. Profiles for a typical steady state banded solution. 
Shown are the local constitutive relation and banded strain 
rate profile, analogous to those in Fig. ^, as well as the com- 
ponents S, W, Z of the polymer stress, as a function of the 
spatial coordinate x. 

We have repeated the four histories with 2? 7^ 0. In this 
case banded profiles are smooth, as shown in Fig. 0. The 
local strain rate j{r) follows that of either the high or low 
strain rate constitutive branches except for the interface 

region whose width scales like \/^. The flow curves are 
shown in Fig. as can be seen, the fiow curves are essen- 
tially history- independent. Our main result is that even 
a small diffusion term lifts the continuous degeneracy of 
steady states. We shall see that this occurs by select- 
ing the value of the total shear stress at the position of 



The history independence and uniqueness of banded 
steady flow is shown by the coincidence of the segment 
Gi — G2 in Fig. |^ for all scenarios and high curvature 
p — 0.1. However, some history dependence remains for 
small curvatures, as seen by the two plateaus G1—O—G3 
and G5-O-G4 for p = 0.01 (Fig. |b). The positive 
slope plateau G1—O—G3 corresponds to the conventional 
two-band flow with the high shear rate band at the inner 
cylinder, while for the negative slope plateau G5—O—G4 
the band sequence is inverted. The normal band sequence 
can be obtained by jumping up from rest for small gap 
velocities (between Gi and O) and the inverted sequence 
can be obtained by jumping up from rest only for high 
gap velocities (between O and G4). In order to scan the 
entire lengths G1—G3 and G4—G5 of the plateaus, ramp- 
ing is necessary after jumping up, as shown in Fig. ||. 
We have checked that both band sequences are linearly 
stable by considering small perturbations of the inter- 
face position and shear rate. Neverthess, as discussed in 
(Radulescu and Olmsted, 1999), the inverted band se- 
quence is only metastable and nucleation processes may 
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change the band sequence. An inverted band is possi- 
ble when both the high and low strain rate branches of 
the constitutive curve span the entire gap, as in Fig. ||c. 
There has been no experimental evidence yet for the exis- 



tence of the metastable inverted band sequence, but this 
could be because either most experiments are performed 
by jumping up at low V, or apparatus noise and interface 
oscillations destabilize the inverted band sequence. 




4.2. Selection mechanism 




0.480 



0.482 0.483 
Shear stress 



0.495 



FIG. 9. Several interface profiles in the representation 7(0-). (a) two different sequences of the bands (high-low and low-high), 
p = 0.01, V = 1.5, 2.5, 3.5, 4.75 and V = 0.001, 0.003; (b) p = 0.1 and f) = 0.001, 0.003, 0.005 and V = 2.0, 3.0, 4.0, 5.0, 6.0. For 
a given "D the profiles for different V differ imperceptibly. 



The key feature of stress selection by diffusion is that 
a stationary interface is stable only at the unique posi- 
tion r» inside the gap corresponding to a selected value 
(7* (as long as is a distance from the walls which is 
large compared to the interface width). A rigorous proof 
of this property does not exist in the case of J-S model, 



but justifications can be found in Lu et al. (199£), and 
Radulescu and Olmstec^ (199!;). To check this we have 
represented the steady shear rate profiles as a function of 
a = T/f^ instead of f (Fig. ^ . In this representation dif- 
ferent banded profiles for different values of V, and hence 
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f*, practically coincide for the same value of V. Profiles 
obtained for different values of V or with different band 
sequences intersect at the same stress &* = a^^i within 
0.01%, which we define as the selected stress (see Fig. ||). 
This selected stress is independent of the cell curvature, 
the diffusion coefficient 2?, and the gap velocity V . For 
example, for e = 0.05, a — 0.5 we find a selected stress 
a*^i ~ 0.4828, which is almost identical to the selected 
stress for the plan ar JS model for the same viscosity ratio. 



0.4827 (Lu_etaL, 1999) 



Hence, for interface widths much smaller than the radii 
of curvature or cylinder gap, the flow is essentially homo- 
geneous on the scale of the interface, and the interface 
migrates to that position in the cylinder for which the 
selected stress cr*g; obtains. The selection of a*^i pro- 
vides a relation F/fJ = (7*^^ = const, between the torque 
and the interface position, thus fixing F for a given gap 
velocity V and band sequence. The magnitude of T) has 
negligible effect in the thin interface limit. 

Stress selection has interesting consequences for the ki- 
netics, one of which is illustrated in Fig. |l^. For flows 
jumping up from rest with gap velocities in the banding 
regime, the kinetics until an interface is well-developed 



are the same as the kinetics for I? = 0, as see n by co- 
incidence of the solid and broken curves in Fig. 10, For 
P = the kinetics ends with the formation of one or sev- 
eral interfaces, while in the presence of diffusion a second 
stage follows in which the interface migrates to its equi- 
librium position at the selected stress. For cylindrical 
Couette flow the stress is a monotonic function of posi- 
tion, so this equilibrium position is unique and steady 
banded flow has only one interface. If multiple interfaces 
form the "excess" interfaces will eventually be expelled 
at the walls. This is the case for low curvature (Fig. foh. 
The three band solution obtained when 2? = (FigT^) 
becomes a transient for "D ^ 0. For low gap velocities 
the interface closest to the inner cylinder is expelled and 
the normal sequence of bands is obtained at stationarity. 
Note that the expulsion time increases for decreasing 2? 
(compare V = 0.0032 and V = 0.005). The inverted 
sequence results from an initial three-band transient at 
higher velocities when the interface closest to the outer 
cylinder is eliminated. The timescale of the interface dis- 
placement is long for small V, so a three-band transient 
could have a very long life and be mistaken for steady 
flow. 
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FIG. 10. (left) Strain rate vs. time for p = 0.01,2? = 0.0,0.0032, and 0.005 and V = 3.0. (right) Torque vs. position for 
T) — 0.005. Elimination of the interface occurs at a local minimum (in time) of the torque. Mesh of 200 space points, time step 
0.00024, for startup from rest. 



4.3. Hysteresis 



Ramping up and down at the extremities of the 
plateaus define different paths, as shown in Fig. the 
low strain rate branch (Fig. pT^,c) or high strain rate 
branch (Fig. [Tl|b,d) and the banded plateaus both pos- 
sess locally stable steady state flows for a narrow range of 
strain rates. Upon ramping up from zero the high shear 
rate band begins to form at a gap velocity Vn (point B2 
in Fig. Hi) which is slightly different than the top-jump 



value. Vn converges to the top-jump value when P — > 0. 
On ramping back down the interface between bands is 
eliminated when it reaches a distance from the wall of 
order its width (point Bi). The difference between the 
corresponding gap velocity 14 and the extrapolation to 
the low strain rate branch scales like VP- Furthermore 
the width of the hysteresis loops, \Vn — Ve\, increases 
with decreasing V, consistent with the banded plateau 
remaining until the interface "touches" a cylinder wall. 
For large enough V hysteresis vanishes. 
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For small V the hysteresis loop at the low shear rate 
end of the plateau contains the metastable part Bi — B2 
of the low shear rate branch. This is reminiscent of the 



experimental results of Grand ct al, ( 1997 ): under con- 
trolled strain rate conditions they found a stress (analo- 
gous to El in Fig. [Tl|a) below which the system remained 



homogeneous, and above which the high strain rate band 
nucleated. Similar phenomena are expected at the high 
shear rate end of the plateau, but an experimental study 
of this region raises difficult problems because of various 
noise sources. 
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FIG. 11. Small velocity step ramping up and down [dV = 0.002) in the hysteresis regions, with diffusion, a) p — 0.1 left end 
of plateau; h) p = 0.1 right end of plateau; c) p = 0.01, left end of plateau, d) p = 0.01, right end of plateau. The spatial mesh 
is 200 points, with a time step of 0.005. In all cases the diffusionless top and buttom jump are shown, respectively. 



5. SUMMARY 



We have studied shear banding in the local {T> — 0) 
Johnson-Segalman model in cylindrical Coucttc flow un- 
der controlled strain rate conditions. The inhomogeneous 
flow conditions of the Couette geometry generate history 
dependent steady states with a single interface of zero 
width. We have outlined some of the behavior under var- 
ious flow histories, and identified families of attractors 
which govern the behavior in the absence of non-local 
terms in the dynamics. 



To recover a unique stress and the history indepen- 
dence seen in experiments, it is necessary to augment 
local constitutive equations with non-local information. 
We have added a simple diffusion term and shown that 
the banding profile is, in the main, history independent. 
The steady state flow curves we have found are consis- 
tent with the growing body of work on wormlike mi- 
celles; moreover, regions of the flow curves suggest the 
analogs of metastability, supercooling, and nucleation 
and growth, and are suggestive of recen t exper i ments 
on wormlike micelles ( [Berret et al , 1994; Berret, 1997 
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Grand et al., 1997). The flow curves have a selected 



banding torque which increases with increasing gap ve- 
locity, with a more pronounced increase for more highly 
curved Couette geometries. 

The diffusion term represents a singular perturbation 
which allows the interface to only lie at that position 
for which the shear stress has a well defined selected 
value. In the Couette geometry this position in unique, 
and multi-interface steady state solutions are prohibited, 
n For flatter geometries a metastable inverted band se- 
quence exists, with the high strain rate phase near the 
outer cylinder. However, there is as yet no accepted the- 
ory for diffusion terms in wormlike micelle solutions, and 
nor have they been measured. Our results suggest further 
interesting effects to study, such as band nucleation and 
expulsion, interface motion, and kinetics, which we hope 
to address in the future. For example, V controls inter- 
face speed, which thus provides a possible experimental 
measure of this unknown quantity. 

Note Added — After this manu s cript was completed we 
learned of recent work by |Yuan| ( |1999| ), in which the JS 
model was supplemented with a diffusion term in the 
strain rate (effectively a non-local viscosity) , and calcula- 
tions were performed using a 2-l-lD Eulerian-Lagrangian 
technique in planar shear flow. Yuan also found unique 
stress selection in the presence of a diffusion term, but the 
results with no diffusive term still gave unique (and dif- 
ferent) stress selection. This indicates that either the nu- 
merical discretization technique introduces spurious non- 
local terms into the dynamics, or that the larger fluctua- 
tion phase space of higher dimensions somehow stabilizes 
a unique stress in the absence of diffusion terms. The lat- 



ter possibility appears to contradict Rcnardy (1995). 
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